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Abstract 
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strong alternating field of high frequency oj is studied in terms of the 
Schwinger-Keldysh nonequilibrium Green functions. The nonlinear 
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resonance structure is elucidated. In particular, it is shown that under 
certain reasonable conditions on the physical parameters, the Fourier 
component /„ is sharply peaked at n = where /S.E is the spacing 
between two levels. This frequency multiplication results from the 
highly nonlinear process of n photon absorption (or emission) by the 
tunneling system. It is also conjectured that this effect (which so 
far is studied mainly in the context of nonlinear optics) might be 
experimentally feasible. 
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1 Introduction 



The physics of resonance tunnehng through a single quantum well is at the 
center of theoretical and experimental activity for more than two decades. 
Electron (or hole) confinement between penetrable barriers in a semicon- 
ductor enables the investigation of numerous important phenomena such as 
negative differential resistance, Coulomb blockade, single electron tunneling, 
single electron pump and many others. The pertinent physics is very attrac- 
tive both because of its richness and its potential device oriented nature. 

At the core of the phenomena of resonance tunneling lies the relatively 
simple picture of a small cavity connected by tunneling barriers to two parti- 
cle reservoirs (also termed as leads). The dimensions of the cavity are small 
enough so that the energy levels inside it are well separated. These levels 
might be either single particle levels whose spacing is determined solely by 
geometrical considerations or a few particle levels determined by interactions 
(as in the Coulomb blockade systems). If there is a difference between the 
chemical potential of the left lead {(II) and that of the right lead, (/i/,) a tun- 
neling current occurs between the two leads. If an energy level of electrons 
in the cavity occurs between fi^ and fiji this current displays a resonance 
structure. Evidently, the pertinent physics is time independent, namely, one 
speaks here of a DC current. Moreover, in most cases, resonance tunneling 
through a single level can be treated within the formalism of linear response. 

Recently, interest is directed toward non linear time dependent transport 
phenomena in double barrier resonance tunneling systems. Experimentally, 
investigation of AC current in mesoscopic devices proves to be feasible [0]. 
(Yet, theoretical analysis of the above experiment remains in the realm of 
linear response 0, |^, ^). 

The relevant physics inevitably becomes richer and more difficult to an- 
alyze. It touches upon qualitatively new phenomena which depend on how 
space and time dependent electronic states interfere. Among the effects 
which are inherently based on non-linear response one might consider elec- 
tron pumps p, H, 1^ , lasers , photon assisted tunneling p, |10|, |ll], |12|, |13|, |14 
and others. 

In the present work we concentrate on a relatively new effect, namely, 
frequency multiplication of current response. As a theoretical model one 
may consider a double barrier resonance tunneling system containing at least 
two quantum levels whose spacing A introduces a new energy scale into 
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the problem. The system is then subject to a strong monochromatic AC 
voltage of strength W and frequency u, which results in a non- linear current 
response I{t), the main object of our study. The combination of strong AC 
electric field and level interaction with the electrons in the reservoirs might 
lead to transitions between these levels which are assisted by many photon 
emission (or absorption), a highly non-linear effect. It can be analyzed in 
terms of the Fourier components /„ of the current I{t). In linear response, 
one expects lo and Ii to be the only non-zero components. Here, however, 
higher components may be significant. As we show below, /„ as a function 
of n is peaked at n = W/A, namely, there is a resonance when n photons are 
absorbed or emitted following transition between the two levels. This kind 
of frequency multiplication is familiar in non- linear optics but, to our best 
knowledge has not yet been investigated in micro- electronics. We give below 
some realistic estimates of this effect. 

As for the theoretical treatment of the above model, we start from the 
familiar tunneling Hamiltonian and employ the Schwinger-Keldysh nonequi- 
librium Green functions, through which the current is calculated in a straight- 
forward way. In order to avoid complications while stressing the important 
role of two level quantum wells we restrict ourselves in this work to non- 
interacting particles. Principally, interaction effects can partially be included 
within screening or Hartrcc-Fock schemes. As we explain in the next section, 
the case of a Coulomb blockade system in which the interaction is reflected 
through the charging energy is qualitatively included in the framework of our 
approach. The physics of non-linear AC response of a general interacting sys- 
tem in a resonance tunneling device contains new effects which go beyond 
the scope of the present work. 

In the following section the problem is formulated, and the various parts 
of the pertinent Hamiltonians are defined and justified. Then, in section 3 
method of solution in terms of the Keldysh Green function is introduced. In 
particular, the free particle Green functions arc written down and the Dyson 
equation for the "lesser" Green function is derived. Section 4 is devoted to 
the discussion of the tunneling current. An exact expression for the current 
is derived in terms of Keldysh Green functions and the origin of resonances in 
the non-linear response is explicitly elucidated. Analysis of higher harmonic 
generation is carried out in section 5, where numerical results are presented 
and the conditions for obtaining peaks at higher harmonics are discussed. 
The paper is then concluded with a short summary. In Appendix A there 
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is a proof that a pure time dependent potential has no physically observable 
effect, while some technical point related to the derivation of Dyson equation 
in the Keldysh formalism is explained in Appendix B. 



2 Formulation of the problem 



We consider a structure where a charge carrier has two barriers in its way 
like the one drawn in fig. 0. The structure is analogous to the Fabry - 
Perot resonator - the motion of a carrier is almost quantized in the central 
region but it still can escape into the leads. Thus the energy levels of the 
central region provide resonances for the transmission of the charge carries 
from one lead to another. Hereafter we will speak about electrons. The 
same type of structure can be made for holes but degeneracy of the valence 
band can complicate the calculation. An arbitrary combination of dc- and 
ac - potentials is applied to the structure (the potential differences are still 
required to be small compared with the Fermi energies). Our main goal is to 
calculate the time - dependent currents in the system. 

Before doing this let us mention possible experimental realizations. Prac- 
tically the appropriate structure can be obtained in several ways: 

1. By putting several layers of different semiconductors having matching 



lattices and different gap widths one on top of another |T^, |T6|, [17 



The electrons then move in the 
direction perpendicular to the layers. As an example, the profile of the 
conduction band for a particular case of a GaAs - Al^ Gai_xAs structure 
is shown in fig. |I|. 



By application of an external gate potential []T0|, |TT], ^ 
0, |36|, ^ 0. The electrons then move along the layers but 

their movement is restricted by the external repulsive potential. The 
gate is shaped such that the electrons move in one direction and two 
barriers are formed in their way. 



3. By placing an impurity in a tunnel barrier 



In order to describe the dynamics of the system we use the tunneling 

The total Hamiltonian is: 



Hamiltonian method 



H{t)=Hf + Hac{t)+HT 



(1) 
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where 

k n p 

is the Hamiltonian of free particles with neither ac-field nor couphng between 
the leads and the central region. Here e^, e„ and ep are free-particle energies 
in the left lead, the central region and the right lead respectively. The energy 
is taken with respect to the bottom of the conduction band in the central 
region. The operators a^, a^, b^, bp, c+, c„ are creation and annihilation 
operators in the leads and in the central region. Furthermore, k (p) are 
momenta in the direction perpendicular to the layers, n numbers the levels 
in the central region (hereafter k (p) refer to the left (right) lead, n, m, n', 
m' - to the central region), summation over n in (^) expresses the presence 
of more than one energy level. The third term, 

Ht = Y1 {TnkCUk + T^C^Cn) + E (^pn&p + (3) 
k,n P,n 

is the part responsible for the tunneling through the barriers, T^j^^J^ are trans- 
fer matrix elements between the leads and the central region. Finally, the 
time dependent part is 

Hac{t) = Wl sin(cut) ■ E «fc «fc + Wr sin(cut) ■ ^ K^p (4) 

k P 

where sin(a;t) are potential shifts of the leads with respect to the cen- 

tral region caused by an external ac-potential. dc - potential shifts are in- 
cluded into the energies ek{n,p)- The arrangement in which Wl and Wr have 
opposite signs corresponds to an application of an ac-bias as it was done in 
the experiments [1^, |2^, |3^ . The choice Wl = Wr describes an application 
of an ac-voltage to the gate electrode superimposed on the central region (see 
the experiments ||T0| , |11|). We notice here that shifting both leads together 
(having the potential of the central region fixed) is equivalent to shifting the 
central region (having the potentials of the leads fixed) since an application 
of a uniform (time - dependent) potential is not observable even if it is ar- 
bitrary strong and arbitrary fast (see Appendix A). The situation Wl = 0, 
7^ corresponds to an application of an ac - voltage only to one barrier 
as it was done in the experiment . 

Our choice of Hac is based on the following assumptions: 



5 



The electrons in the leads respond to an applied field very fast since 
we deal with the frequencies much less than the plasma frequency. It 
means that any change of the external potential causes an immediate 
rearrangement of the electrons. In other words, the internal potential 
responds very quickly to an external field 
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The concentration of the electrons in the leads is high enough to screen 
an external field |]45| , Therefore potential is uniform in the leads 



and drops in the barriers |, ||, ||, g^, |8|, |9l 0, |T|, |53|, |5|, |5[ 



3. We used a widespread assumption |^ ^ Q that the probability 
of direct transitions between the energy levels in the central region due 
to the ac-field is small and can be neglected. 

We do not restrict our consideration to the case of small Wl, Wr (linear 
response). They can be arbitrary large. A strong ac-field {Wl ^ or 
Wr S> w) leads to a non - linearity. One of its signatures is the generation 
of current harmonics with frequencies much larger than uj. 

The Hamiltonian (|I]) does not include Coulomb interaction. It is fre- 
quently omitted in the treatment of quantum wells ^ |SD|, ^ 



6^ , |63| , |64| , and others] . The number of electrons there is very 
large so the Coulomb interaction can be treated in the mean - field approx- 
imation. In fact it makes concentration - dependent corrections to the one 
- particle energies e„. It can lead to some changes in the dependence of the 
current on the bias (such as the appearance of a hysteresis ||16], ^ |67|1 ) 
or, may be, on the magnitude of the ac-field but the physics of the elec- 
tron - photon interaction remains essentially within the independent particle 
model. In quantum dots the energy of the Coulomb interaction is much 
larger than the spacing between one - particle energy levels, level widths and 
temperature. The dc-bias, the magnitude of the ac-potentials {Wl, Wr) as 
well as the frequency are in general much smaller than the Coulomb interac- 
tion energy. Hence, to a good approximation, every tunneling event involves 
only that electron which happens to have the highest energy. The Coulomb 
repulsion with the others only provides it with the potential energy needed 
to penetrate into the leads. Thus they serve as a background. We can forget 
(at least, qualitatively) about this background in our consideration of the 
interaction of that electron which provides the current with the ac-field. 
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Rigorously speaking, the Hamiltonian (|1|) is one - dimensional. In QWs 
there are lateral degrees of freedom provided by the motion along the layers. 
But at every particular value of the lateral momentum the problem is one - 
dimensional. Even if the dependence of electron energy on lateral momentum 
is important the total current can found by simple integration over it. The 
lateral degrees of freedom can be removed by a magnetic field perpendicular 
to the layers. 



3 Solution in terms of Keldysh Green's func- 
tions. 

In order to work out the current we first have to analyze the electron propa- 
gation. To this end we employ the non-equilibrium Green function technique 
suggested by Schwinger ||68|, KadanofF and Baym |^ and Keldysh fT^ (for 
a review see Ref . []52|, [7l| , ff^ , |T3| , ^]). Since no equilibrium is required one 
can rigorously consider large perturbations and high frequencies drawing the 
system far from its steady state. The method uses a time variable defined on 
two sides of the real axis. It is equivalent to introduction two independent 
Green's functions, one which characterizes the dynamical properties of the 
particles the other one describes their distribution . 

It is convenient to use the following two Green functions (the others can 
be expressed through them): 



Gljih, t) = -lOiti - t) (v9.(ti)<^](t) + 



(5) 



where (pi{ti) (^(p^j{t)^ are operators in the Heisenberg picture representing 

o-kiti), Cn(ti), bp{ti) (cifc(t), c|j(t), &p(t)) in correspondence with what val- 
ues the indexes i and j take {k, n or p). The first one is the usual retarded 

Green function. The advanced Green function (G'^j^ is its Hermitian conju- 

r 1 * 

gate that is: Gfj(ti,t) = G^j(t,ti) . The other one contains information 
about the distribution of electrons and their correlations. At ti = t it gives 
the one - particle density matrix (p) |7T|: Gfj(t, t) = iNpij(t), where is the 
number of particles in the system. The advantage of the method is that it 
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treats G"" and in a unified manner. Both of tliem can be found through 



the Dyson's equation 73 



G"' = /[l + S^'GH (6) 
G< = [1 + G'Sq ^< [1 + + G"S<G'^ (7) 

Here multiphcation imphes summation (or integration) over space variables 
and integration over time, and are Green functions defined as (^ but 
for an unperturbed Hamiltonian, S^, and are proper irreducible self - 
energies. 

An application of the formalism to tunneling systems is especially power- 
ful if one chooses Ht as a perturbation [^). Then then self-energies assume 
a very simple form: 

= Ki^' = 



yr{a) 
^kk' 


- ^pp' 


yr{a) 

— ^kp 

rpL 

— -'-nk 


v-ir 
^nk 
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^kn 
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rpL 

— ^kn 
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pn 


— ^np 

— ^pn 

= 


_ 

np 

_ 

pn 



^nk 



np 



(8) 



This is due to the fact that any vertex caused by Ht has only two entries. 
The only diagram contributing to an irreducible self - energy is the simplest 
one drawn in fig. |^a. The next one drawn in fig. ^ is already reducible. 

In the next section we will show that tunneling currents can be expressed 
solely in terms of the Green functions in the central region. To find them we 
iterate (^ and get: 



Glmiti^t) = g:^{tut)+ (9) 

/ + 00 i" + 00 
/ dt2dt3 gnnij^l, t2)Xnn'{t2, i^3)G'J^'m(i^3, t) 
-oo J — oo 



Xnn'ih,^) = X^^,{t2,t3) + Xj^^,{t2,t3) (10) 

^nn'it2,t3) = "^T^kdlki^^, hWk^' 
k 

^nn'(^2,^3) = ^T^pgppih, t3)T^' 
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Equation establish an algorithm for calculating G^^ without an 

iteration but nevertheless it contains eleven terms on the right hand side. We 
showed (see Appendix B) that only two of them remain after a long enough 
period of time has passed from the moment the tunneling was switched on: 

/+00 r+co 
/ dhdh G:At„t,)Y^,^,{t„ h)Gl,^{t„ t), (11) 
-oo J — oo 

n'm 

where 



k 



V 

We notice that does not depend on the initial distribution of electrons 
in the central region (because g'^^ do not appear in formula (0)). It depends 
only on their distributions in the leads through g'^y. and (7^. 

It is easy to show (see Appendix A) that even under a time-dependent 
potential of an arbitrary large amplitude and an arbitrary high frequency the 
Green functions for electrons in isolated leads are given by "adiabatic - like" 
expressions: 

9'kk\t2M) = -^S,,.e{t2-h)e-'^^s^''^'-''-^'^^^^''= (13) 

9^k'ih,h) = ^5,../L(e,)e^^^'3[^^^^-^*)]-'* = 

= z5fefc'/L(efc)e-^'^(*2-*3)+*^('°''('^*2)-cosM3)) 



provided that the time - dependent perturbing potential (ifac) is uniform in 
every lead. Here [. . .]kk denotes matrix element, fii^-k) = e^-^j^ — is the 

e fee +1 

Fermi function for the left lead, /i^ is its chemical potential, G is the tem- 



perature. Formulas for g^^,, g^^, can be obtained from ([TsD by replacements 
k, k' p,p', Wl Wr, fiL ftR. 

Now we substitute (|13|) into ([lOD and ([l^) and replace the sum over k, p by 



an integral over energy. Further, we employ the "wide - band approximation" 



[ [T^ , [78[] assuming the elastic coupling to the leads 



ri, = 2'np{e,)T^,Ttr,. and F^, = 21, p{e,)T^T^,, (14) 
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to be energy independent (i.e. they do not depend on A;, ek) within the band. 
The band widths and the Fermi energies are large compare to Wl, Wr. Thus 
(|TUp and ([T^) take a rather simple form: 



Xnn'{tl,t2) = --r„„/(5(ti - t2) (15) 

JVl ^tt 

JVr 2,71 

Hereafter Vl (Vr) are energies of the conduction band bottoms in the left 
(right) lead, r„„/ = F^^, + r^„,. We can not put Vl, Vr —oo in the formula 
(|16D because it would lead to a divergence in the ac-current. 

We want to emphasize that substituting (^) into @ yields a time - 
independent equation. The problem then reduces to the task of finding the 
usual retarded Green function for a time - independent potential well with 
finite barriers. This is due to the energy independence of F^^f ■* . The Fourier 
transform of the equation is: 

1 

n' 

A rigorous solution for a two - level system shows that the mixing terms 
^nm)''^ 7^ gi'^^ ^ contribution to the current of the order of F/ |e„ — €^1- 
We assume F ^ l^n — Cml and drop them out hereafter. With the same 
accuracy are given by: 

GU^) = ^— r7 (18) 

where F„ = F„„ + F™, F™ are intrinsic level widths due to inelastic interac- 
tions, leakage of electrons into lateral directions etc. This is equivalent to the 
widely - used assumption that the quasilevels of a potential well with finite 
barriers have complex energies e„ — i^. We see that the main contribution 
of tunneling in the dynamical properties is level broadening. 
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Substituting ( ]T6D into (|TT]) gives an exact expression for G"^ as soon as 
we know C: 



-^^ \ ^—isujti+iqujt 



s,q= 



+ 



Vr 27r 



G;„,(efc + s^)G:;:„,(efc + ga;)+ (19) 



s,q= 



UJ 



UJ 



We notice that the transition amphtudes G'^^.n ^ m are important as we 
show in the next section. 

This difference between GJ^„ and G^^, n ^ m arises from the difference 
between X and Y in the formulas (|15]), (plGl). The quantity X has such a 
simple form because f{ek{p)) does not appear in g^i^, g^^ so the summation 
over k{p) (equivalent to integration over efc(ep)) in formula (p!0D gives a S - 
function ^(^1—^2) if the band (i.e the range of ek{p)) is wide. The occurrence of 
f{ek{p)) in g^i, (g^p) leads to an upper limit of integration in (|1^) around ^l{b)- 
If fiL{R) is large enough the 6(ti — ^2) is recovered. Then the whole solution 
becomes invariant under time translation and the transition amplitudes G^^, 
n ^ m are small, of the order of r/|en — Cml- 



4 The current 

A double barrier structure is integrated within a circuit. The measured 
current is determined by its influence on the circuit. If the barriers are 



modeled as capacitive - resistant elements |^ the currents in the leads are 
given by |T^: 



hit) = ^^^^CLUjWLCOs{uJt)-^CRUjWRCO^{ujt)+ (20) 

where II is the total current in the left lead (t) and /^(t) are tunneling 
currents through the left and right barriers respectively, Cl{r) and Cg are 
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capacitances of the central region relative to the left (right) lead and a gate 
electrode (or another background), C = Cl + Cr + Cg. The current in the 
right lead is obtained by replacing L by R. The first two terms describe the 
contribution of capacitive currents (due to the presence of accumulation and 
depletion layers). They have the frequency of the input ac - voltage. But the 
tunneling currents contain also higher harmonics as we show below. 

The tunneling current from the left lead into the central region is defined 
by the change in the number of electrons in that lead: 



/ dNj 
\ dt 



rpL * 



■ mk 



{alcm) 



X/fcm I^mk^km^i ^) 



(21) 



where Nl = Hk O't^k- The tunneling current from the right lead into the 
central region (/^) is obtained from the above expression by a change L ^ R, 
k ^ p, ak ^ bp. 

Using Dyson's equation for G^^ allows one to express the tunneling cur- 
rents through Green functions in the central region (which we found in the 
previous section): 



(22) 

Under the assumption of energy independent coupling to the leads which 
we used in the previous section (see (0)) we can use the formulas (|T5|), ([16|) , 



(|T8|) and (p!9|). Dropping out the terms that are always small of the order of 
r/ \ej — tj'l in comparison with the others we obtain after some algebra: 



A + E \Il\ cos(ri^t + 02), where02 = arg 



ra=l 



2e 



- ^ [A,{n) + A]{-n) + i?,(n) + B*{-n) 



^ n' 



nw + 1 



+ 0O 



3 I " 2 



(23) 



E Js{—)JsU—)F 



ijj 



UJ 



5 ■> 

IxJ UJ 
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J' €j' — Sj — nuj + i 



ri+r, 




'e,- -11 _ ii-V ^ke\_l_ f+oo f{e)de _ 
cu CO CO 2Tr Jv f -l su: - f , + r-^ 



\ oj UJ UJ UJ ) 2t: Jv e + suj-ej + i^ + 1 

Analogous expressions hold for /^(t) as well. For the rest of our discussion 
we choose to study /J. We point out that the current consists of many 
harmonics with the frequencies nu. The sums over express summation 
over all energy levels in the central region. We notice that the terms Aj{n) 
and Bj{n) are quite similar in structure. The first one includes the currents 
due to direct photon - assisted transitions between the left lead and the 
central region and due to the effect of the left lead on the population of the 
central region. The second one describes the current through the left barrier 
due to the influence of the right lead on the electrons in the central region. 
The main difference between them (beside a few replacements of L into R) is 
that the first term in the brackets in A (i.e. i^^j) is absent in Bj{n). At zero 
temperature the integral in the expression for F can be easily computed: 

ej- li _ ^ H-V Vj ^\ _ I -suj-ii-i^ 



s,- ,^,0 =— In 



UJ UJ UJ J 2tt ej - suj -V - i ^ 

The rest of our calculation is done at zero temperature. 

The term in A containing Tjj gives an additive contribution of different 
energy levels to the current. The same is true for the terms of the sums over 
j' (in Aj{n) and Bj{n)) with j' = j. The other terms in these sums are caused 
by the presence of different energy levels together. Every one is caused by a 
certain pair of energy levels. It is obvious that they have a resonance nature. 
The j'-th term can be large when |e^/ — Cj — nuj\ < (Tj/ + Tj)/2. Thus it can 
give large ac-current at n^cs ~ ^ ^jl (under condition W > {ej/ — ej\). 
If UJ <ti \ej' — ej\ one has n^es ^ 1, i-e. the frequency in the output is much 
larger than in the input. We notice that this term generates the rires-th 
harmonic exclusively. The others are small of the order of „ , (we 

suppose high frequency F/uj <^ 1). 
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The resonance is caused by many - photon absorption due to the strength 
(W ^ uj) of the ac-field. The formula (|2l|) expresses the current at the 
moment of time t through the probabihty to create an electron in the state 
m annihilating it in the state k as it is drawn in fig. |^a. To do that it is 
necessary to have an electron in the state k and a hole in the state m. The 
formula (^) shows that it is provided by creation of an electron in the state 
k and a hole in the state n at an earlier moment of time ti, ti <t (see fig. Qd). 
In the presence of a strong ac - field the process in which n ^ m should be 
considered. The hole then has to propagate from n to m. The formula (pj]) 
expresses the transition from n to m through the kind of process drawn in 
fig. ^ for ^3 < ^2- An electron is transfered at the moment t^, <t from the 
state m to the state p in a lead leaving a hole in the state m. At the moment 
^2,^2 < ti the electron tunnels to the state n. Eventually an electron is found 
in the state n at ti and a hole is found in the state m at t. In the absence of a 
time - dependent potential this process is impossible because the transitions 
from m to p and from p to n have to conserve energy (after integration over 
t2, ts). Thus m has to be equal to n. Hadf) leads to a multiplication the 

wave function of the state p by the phase factor e'^~'^°^^ ' (see Appendix A 
and formula (|13|) in the main text). In the energy representation the wave 
function becomes dispersed over all energies ep±suj, s integer with the weight 
Js(^). In other words the spectral function of the state p has extra peaks 
with energies ep±suj. Their magnitudes are proportional to (^) (see 0). 
Thus it is possible to have transitions from m to p and from p to n when 

Cm - ep = gtu, e„ - ep = suj (24) 

This is the origin of Jg and Jq in the formulas (|19|) , (|23D . These transitions are 
accompanied by emission (absorption) of q (s) photons. The total number 
of absorbed photons is s — q = n. When ac-field is strong (Wr ^ u ) 
Jsi—) and Jg(— ) can be large for s,g ^ 1. This leads to a possibility of 
many - photon absorption (emission). The resonance conditions (^) give 
rires = s — q = '''"J'^" . It is easy to show that, in fact, Wr > Tires'^ is required 
to obtain large nres-th harmonic. The physical content of fig. |^c expresses 
the contribution of ^^(^2,^3) into (|lID . The same process but through the 
left lead expresses the contribution of F^(t2,^3)- At ^2 < the picture is 
very similar. The process just starts from a transfer of an electron from p to 
n at ^2- Then the hole tunnels from p to m at ^3. 
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If we do not assume that the lead energy bands are wide (in particular 



that 



L(R) 



is large), then the only change at Wl = (which is required in 
order to observe the high frequency resonance as we show in the next section) 
is that now G^^^^n 7^ m is not negligible. These non-diagonal elements of 
the retarded Green functions could make a positive or negative contribution 
to the process drawn in fig. ^jc. Yet, we argue that the resonance would 
persist also in that case. 

Every pair of energy levels in the central region gives an independent 
contribution to the current (^). In the next section we consider one pair of 
energy levels and explore the dependence of the current it produces on the 
parameters of the system. 



5 High harmonics generation. 

In this section we consider a strong ac - field {Wl ^ u; or Wr ^ u) and 
show the dependence of the current on the parameters of the system. Since 
every pair of energy levels gives a separate contribution to the current it is 
enough to consider only one pair. 

Generally a strong ac - field leads to generation of high harmonics. In 
fig.^we plot the amplitude of the harmonics (|/£|) versus their number n for 
(as an example) Wl = Wr = 30, ei = 11.5, €2 = 30.5, fiR = 20, /i^ = 5, 
Vr = Vl = -50, T^l^^ = Tai^^ = Tf^ = -Vf^ = 0.05, 00 = 1. We did not 
show the dc - component because its magnitude is much larger. We did not 
show the first harmonic (whose amplitude is also a few times larger than that 
of the others) since it is anyhow only a part of the current response with the 
frequency u. The other part is given by the capacitive currents (first terms in 



the formula (|20|) . The main contribution to the broadened spectrum plotted 
in fig.^ is due to the direct photon - assisted transitions from the left lead 
to the central region and back described by the first term in the braces in 
Aj{n) (see formula (^)). The terms with j 7^ j' arising from the presence of 
two or more levels together can lead to generation of one solitary (and high) 
harmonic (see fig.§). In the rest of this section we describe the conditions 
required to observe one high harmonic alone and show its dependence on the 
parameters of the system. 

Since the first term in Aj{n) generates many harmonics it must be elimi- 
nated. The easiest way to do it is to apply ac - voltage only to the right 
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barrier: Wl = (in fact it is enough to have Wl <^ |e2 — ei|)- Then 
Aj{n) = 0, n^O because J, (^) = J,(0) = 0, s ^ 0. 

In fig.^ we plotted the amplitude of the harmonics (|/£|) versus their 
number n for Wl = 0, Wr = 30, ei = 11.5, 62 = 30.5, /i/j = 20, fiL = 5, 
Vr = Vl = -50, r^^^^ = Tai^^ = = -Ffa = 0.05, to = 1. As before we 
did not plot here the the dc - current and the first harmonic. They are are 
of course large compared with higher harmonics. We notice however that if 
Wr becomes larger than the Fermi energy in the right lead (/iR — Vr) they 
are substantially reduced, and the amplitude of the resonant harmonic might 
be comparable to the magnitude of the dc - current. It is clearly seen in this 
figure that among the high harmonics only the 19-th one {nj-^s = = 19) 
is generated. 

At Wl = and Fi r2 a simpler expression is easily obtained for the 
resonant harmonic: 



4--^^|[5,(nre.) + 5,(n,es)]-r"--|^^ t Js{^)Js.n..A^) 

(25) 

Only the terms with s so that ei — su is inside the right band: Vr < ei — suj < 
fiR contribute to the sum. We notice that Js(^) tends to zero when s 
becomes larger than ^ if ^ 3> 1 so Vr? does not influence the results as 
soon as Vr < ei — Wr. 



In fig.^ we drew the dependence of the resonant harmonic 



r^res 



on 



both ^ (i.e. ac - voltage) and ^J— ^ (determined by the dc - bias) for 



<;2— gi 



8, Fff = Tsf = ^12/^ = -r^i/^ = 0.05, fiR/u = 20, 



Vr/uj = —200. There is no generation (^I^^^^ = 0) if Wr < ni.esC^/2 = 
(e2 — ei)/2 because the conditions |s| < ^ and \s + Uj-esl < ^ can not be 
satisfied together so one of Js(— ), Js+n^^^ (lii.) is small. The generation 
becomes significant for Wr > 62 — ei. It is especially important if there 
are more than two energy levels: distant pairs of levels do not generate 
harmonics. 

The dependence on the bias (i.e. ^^-^f^) is oscillating. The reason is the 
interference of different components of the wave function (remember that in 
a strong ac - field it is spread over a set of energies ± su with different 
s). It shows in the formulas as oscillating behavior of Js(^) via s. Let us 
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consider the sum in (^) at |ei 
to divide it into two parts: 



finl ^ Wr and Vr < ei — Wr. It is useful 



,Wr 



S = Sl 



:Wr. 



UJ 



J, 



UJ ^^^^ UJ UJ 



where Si = J^^ and -C Wr. The second part does not depend on 
ei — fiR. The Bessel's functions in the first part can be approximated by cos: 



J. 



•Wr, 



UJ 




2 I V -n 

cos \ V — s arcsm ; 

-nv \ Wr UJ 4 



since Isl <^ — . Then 



V 



2e 2rf'2r|^]^ 



^ Ti + ttWr/uj 
1 . 



S = Sl 



COS a 



2e 2rf'2r^]^ 



h Ti + r2 Tin,- 



where a 



2Wr 



and [3 



sm a 



ei 



Wr ' 



W. 



+ cos(/9 — ns) + const 



const) (26) 



7r(n,-cs +1) 



. When ei—/iij changes from 
—Wr to Wr the value /^'''^'^ oscillates about rires/Tr times. The amplitude of 



the current 



r^res 



has about 2nres/7r maxima (in fig.|]the number of maxima 
is even larger because the frequency of oscillations is larger at ei — ^r ^ 
±.Wr). Notice that the current is not generated ii \ei — hr\ > Wr. It is 
significant for systems with many energy levels: only those pairs of levels 
that are in the energy range from fiR — Wr to fiR + Wr generate harmonics. 

The dependence of the current on the transparency of the barriers is clear 
from formulas (p3D, (p5|). In fig.^ we showed the spectrum of the tunneling 
current at T/uj = 1/20. If F increases two sets of harmonics grow in the 
vicinity of the 0-th and the rires -th harmonics. The ratio of the amplitudes 
of the side harmonics to the amplitude of the leading (0-th or n^-es -th) one is 
about r/nuj and T/ {n—rires)^ respectively. On the other hand, the magnitude 
of the resonant current is proportional to yI+V2 ' 

The dependence of the current on frequency is oscillatory. If a rectifying 
device like a diode is placed in the output the harmonics contribute to the 
dc - current. In fig.|^ we draw the value (/) = | + ~ Y^n=i I-^lI coming 
from time - averaging of lj^{t) as a function of the input frequency uj at 
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Wl = Q,Wr = 30, ei = 11.5, £2 = 31.5, ^lR = 20, /i^ = 5, Vr = Vl = -200, 
■pL^(i?) _ p^^C-fJ) _ pL^ _ — = 0.05. The peaks correspond to different 
harmonics (n = 22, 21, 20, 19, 18) satisfying the resonant condition n = . 

We have discussed the high harmonic resonance in the tunnehng current 
through the left barrier. To observe it the measuring device must be made 
sensitive to the tunnehng current only through this barrier. Otherwise other 
harmonics will mask the resonance (like in fig.^. If the device measures the 
current in the leads we propose two ways to obtain it based on the formula 

1. Making the gate capacitance small: Cg <^ Cl,Cr and the capacitance 
of the right barrier much larger than the capacitance of the left one: 
Cl ^ Cr. It can be achieved by making the left barrier thicker than 
the right one. Notice that the right barrier must be made slightly 
higher so that and Tr are of the same order (then /J and are 
of the same order). Thus the high harmonic currents in the leads are 
approximately given by: 

hit) - ii{t)-^im^im 

We omitted here the contribution of the capacitive currents having 
frequency uj. The high harmonic currents in both leads are determined 

T 
L ■ 



Making the gate capacitance large: Cg ^ Cl (the ratio of barrier 
capacitances Cl and Cr is arbitrary, there is no need to make one 
of the barriers higher than the other one). Then (dropping out the 
capacitive currents) we obtain from (^0]) : 



hit) - ii{t)-^,^im 

^R + 

^R + + 

We notice that in this case high harmonics in the left lead are produced 
by /J but in the right lead by both /J and (unless Cr 3> Cg). If 
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Wl = 0, Wr 7^ the term A generating a wide spectrum of harmonics 
vanishes in but it is not zero in the expression for /^(t). Then the 
current in the right lead lR{t) consists of many harmonics. To detect 
a single harmonic the current in the left lead should be measured. 

6 Summary 

The basic physical problem addressed in this work is concerned with a non- 
linear response of a two(or more)-level system to a high frequency external 
field. This makes the formalism developed above particularly attractive since 
two-level systems are an important model for many realistic physical situa- 
tions. The most familiar one is evidently a two-level atomic system, whose 
response to a laser field is one of the hall marks of non-linear quantum optics, 
one of whose signatures is higher frequency generation. Here we have focused 
on an electronic analog, where the response is a tunneling current instead of 
an emitted hght. Note that, unlike the optical analog, the theoretical formu- 
lation requires computation of non-equilibrium (Schwinger-Keldysh) Green 
functions. As far as the experimental situation is concerned, we see no real 
obstacle in the road for actual observation of this effect. 
We expect new effects to occur when the many-body physics is included. 
When the size of the dot is small enough then, at voltages and frequencies 
less then the Coulomb interaction energy (which is frequently the case), it 
can be considered as an Anderson impurity, which shows a Kondo type effect 
if the resonance level is deep below the Fermi level and the temperature is low 
enough. Hence, the first problem which comes into mind in this context is a 
non-linear response of a magnetic impurity to a time-dependent field. The 
Kondo effect out of equilibrium has been studied recently by several authors, 
but this is done primarily in the non-crossing approximation, which is valid 
much above the Kondo temperature. The formalism developed above gives 
the hope that also crossed diagrams can be included (the first ones appear 
when a sixth order term in the tunneling matrix elements are computed) 
and therefore it can approach the physics at temperatures lower than the 
Kondo temperature. If external voltages or frequencies arc comparable with 
the Coulomb interaction energy mutual time-dependent resonant tunneling 
of electrons might show new physics. We suppose that this problem can be 
handled using the finite - U Anderson model. 
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Appendix A 

In this Appendix we show that an arbitrary strong time - dependent potential 
has no effect if it is uniform in space. We also calculate here Green functions 
for an isolated lead. 

If the alternating potential {Hac{t)) is uniform the Hamiltonian can be 
written in the following form: 

H{t)=Hf + Hac{t) (27) 

where Hf is time - independent, Hadt) is space - independent. The solutions 
of the Schrodinger's equation are: 

vi>fc(t) = ¥,fc(t)e-^/'^-(*i)'^*^ (28) 

where ipkif) are eigenstates of Hf. The time dependent part Hadt) influences 
only the phases giving all the solutions the same phase factor e~« J Hac(ti)dti^ 
The phase differences are determined by Hf only. Thus Hadt) has no physical 
effect. This simple result is the time - dependent analog of the fact that the 
(time - independent) reference point of energy can be arbitrarily chosen. 

We used this result to write the ac - part of the Hamiltonian of a double 
- barrier structure in the form (^). The ac - shift of the central region is 
ignored. 

In an isolated lead the ac - potential is uniform (see (^). Then the 
evolution operator {U{t^t')) is obtained from the equation ih-^U{t,t') = 
H{t)U{t,t') with H{t) given by (pTj). Its matrix elements are: 

Ukkit,t')=6kk'e-^^^'"'"'^''^'''' (29) 

Skk' appears because H^k' = 0, k ^ k' (i.e a uniform potential does not cause 
transitions). Using the definition of Green functions in the Schrodinger rep- 
resentation with the evolution operator ( pQ]) it is easy to obtain the formulas 
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Appendix B 



In this Appendix we analyze the Dyson's equation for and obtain the 
formula (|Tll) . The Dyson's equation holds for G^^: 

^nm ^nm(.9jj^ ~^ -^nmiydkki 9pp) (^^) 

^nmidjj) = 9nm ~^ ^nk'^kmQmm ~^ ^np'^pmdmm ~^ 

~^ dnn^nkGj^^ + Q^^T^pG^j^ + 

+ ^nk'^ki9u'^ik'G'^,^ + Gl^pTpiQ^TipiGp,^ + 

"I" ^nk'^kiQii'^ipGpra + ^ np'^piduTikG 

^nmidkky 9pp) = ^nri''^n'k9kk'^km'G^(^ + G^^^iTnipQppTp^iGj^i^ 

Here multiplication implies integration over time and summation over re- 
peated indexes; j = m,n,i belong to the central region, k, k' — to the left 
lead, p, p' — to the right lead. 

^nmidfj) depends on the initial state of the central region. It does not 
depend on g^k^Opp {G'', G"^ are determined by the equation (||) which does 
not depend on gkkiQpp)- It describes the behavior of the system with empty 
leads. Indeed, if = ^.Qpp = we have F^^ = 0. Then G<^ = F^^. 
When the leads are empty the central region empties with time. tends 
to zero. After a long enough period of time has passed from the moment 
the tunneling was switched on it can be neglected. Then G^^ = F^^. This 
is the formula ([lT|). We can say that describes the transient processes 
while F^^ gives some kind of a quasistationary (but fully time - dependent) 
solution. 
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Figure 1: Schematic energy diagram of the conduction band for an 
AlxGai^xAs quantum well in the absence of an ac-field. x is the axis perpen- 
dicular to the layers, /xx, and /ir are the left and the right chemical potentials 
correspondently, Vl and Vr are the energies of the conduction band bottoms 
of the leads, ei and €2 are the energies of the levels in the central region. 

Figure 2: Tunneling diagrams. A cross denotes a tunneling event, a thing 
line - free propagator, a thick line - full propagator, a) Dyson's equation. 
Only a single tunneling event contributes to the self - energy, b) any other 
diagram is reducible. 

Figure 3: Tunneling processes leading to the high harmonic resonance. A 
thick black line denotes full electron propagator, a thick white one - full 
hole propagator, a thing black line - free electron propagator, a dashed line 
- tunneling event, a) The tunneling current is produced by electrons going 
from the lead to a level in the central region and back, b) This process can 
be assisted by other levels if ac - field is strong, c) Resonance transfer of an 
electron from the level m to the level n possible in a strong ac-field. 

Figure 4: Spectrum of the tunneling current (amplitude of the harmonics via 
their number). General situation. Wl — Wr — 30, ei = 11.5, 63 = 30.5, 

i^R = 20, ^5,Vr^Vl^ -50, r = rff ) = rg''^ = rf^ = -rf^ = 0.05, 

cu — 1. The dc - current and the 1-st harmonic are not shown. 

Figure 5: Spectrum of the tunneling current (amplitude of the harmonics 
via their number) when the alternating field is applied on the right barrier: 

Wl = 0,Wr = 30, ei = 11.5, £2 = 30.5, fiR = 20, /i^ = 5,Vr^Vl = -50, 
r = T^i^^ = Fsi^^ = rf'2 = -rfa = 0.05, u^l. The dc - current and the 
1-st harmonic are not shown. 

Figure 6: Dependence of the resonant harmonic jj^^res ~j ^^le dc - bias 
(i.e. 61 — ^r) and the amplitude of the ac - voltage Wr. Uj-^s = 8, Wl = 0, 

62 = 61+8, /x« = 20, Vr = -200, r^-^^) = r^f ) = rf^ = -rf^ = 0.05, ^ = 1. 



Figure 7: The output current rectified by a diode versus the frequency on of 
the input ac - voltage. Wl = 0, Wr = 30, ei = 11.5, 62 = 31.5, fiR = 20, 
I^L = 5,Vr = Vl = -200, F = Ff(^) = Fg^) = Ff^ = -Ff^ = 0.05. The 
peaks correspond to different harmonics {n — 22, 21, 20, 19, 18) satisfying the 
resonant condition n — . 28 



